--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\matt\publications\ajps4\replication\data_analysis\individual_religious_participation.log
  log type:  text
 opened on:  11 Oct 2012, 05:39:54

. #delimit ;
delimiter now ;
. *     ***************************************************************** *;
. *     ***************************************************************** *;
. *       File-Name:      individual_religious_participation.do           *;
. *       Date:           October 11, 2012                                *;
. *       Author:         MRG                                             *;
. *       Purpose:        Take individual_religious_participation.dta     *;
. *                       and replicate the results in Table 1            *;
. *           Input File:     individual_religious_participation.dta          *;
. *       Output File:    individual_religious_participation.log          *;
. *       Data Output:    none                                            *;
.              *       Previous file:  individual_religious_participation.dta          *;
. *       Machine:        office computer                                                 *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. set mem 400m;

Current memory allocation

                    current                                 memory usage
    settable          value     description                 (1M = 1024k)
    --------------------------------------------------------------------
    set maxvar         5000     max. variables allowed           1.947M
    set memory          400M    max. data space                400.000M
    set matsize         400     max. RHS vars in models          1.254M
                                                            -----------
                                                               403.201M

. use "C:\matt\publications\ajps4\replication\data_analysis\individual_religious_participation.dta", clear;

. set more off;

. sum;

    Variable |       Obs        Mean    Std. Dev.       Min        Max
-------------+--------------------------------------------------------
        year |    264370    1994.946     5.85372       1981       2004
       ccode |    264370    383.2626    221.2662          2        920
     country |         0
attend_rel~s |    249063    4.335658    2.558971          1          8
        wave |    264370    2.951246    .9936813          1          4
-------------+--------------------------------------------------------
highest_ed~n |    184811    4.433659     2.27252          1          8
education_~e |    184811    1.844944    .7394981          1          3
social_class |    131988    2.667371    .9835625          1          5
income_scale |    225984    4.682708     2.47316          1         10
income_level |    224694    1.954391    .7923673          1          3
-------------+--------------------------------------------------------
        male |    264175    .4814157    .4996555          0          1
         hdi |    257484    .8018018    .1191805       .411      .9486
government~n |    264370    3.873803    2.853166          0   9.444502
social_reg~n |    264370     4.05382    3.079057          0         10
government~m |    264370    2.405628     .829125          1          3
-------------+--------------------------------------------------------
relfrac_al~a |    264370     .439457    .2321796      .0035      .8603
pr~t_alesina |    264370    17.61932    27.12147          0      89.29
roman_cath~a |    264370    33.50454    35.75483          0      95.22
muslim_ale~a |    264370    15.14659    31.12567          0      99.93
taoist_ale~a |    264370    .0533434    .9806157          0      18.08
-------------+--------------------------------------------------------
buddhist_a~a |    264370    1.913725    8.396807          0      66.66
inequality~i |    244467    40.99867    9.505983   17.77497     61.819
urbanization |    261265     65.6716    17.03078         12        100
  regulation |    260589    2.093953    3.418767          0         13
         gir |    260589    23.08274     14.9805          0   77.55682
-------------+--------------------------------------------------------
      region |    264370    7.419272    2.846142          1         10
     polity2 |    239936    6.556536    5.096141         -7         10
     gdp1000 |    215276    12.57108    7.812932    .706843   41.35404
   communist |    264370    .0448954    .2070748          0          1
postcommun~t |    264370    .2380452    .4258877          0          1
-------------+--------------------------------------------------------
pr~y_alesina |    264370    .1369596    .3438053          0          1
catholic_m~a |    264370    .3396338    .4735858          0          1
muslim_maj~a |    264370    .1266331    .3325622          0          1
        lhdi |    257484   -.2341237    .1708629   -.889162  -.0527681
    lgdp1000 |    215276    2.273898    .8071411  -.3469467    3.72217
-------------+--------------------------------------------------------
       age65 |    261392    .1092306    .3119289          0          1
 sub_saharan |    264370    .0637932    .2443846          0          1
  south_asia |    264370    .0465257    .2106211          0          1
   east_asia |    264370    .0513712    .2207541          0          1
     se_asia |    264370    .0223777    .1479089          0          1
-------------+--------------------------------------------------------
     oceania |    264370    .0169346    .1290267          0          1
 middle_east |    264370    .0865113    .2811182          0          1
latin_amer~a |    264370    .1143662    .3182562          0          1
north_amer~a |    264370    .0447138    .2066753          0          1
 east_europe |    264370    .2772743    .4476539          0          1
-------------+--------------------------------------------------------
 west_europe |    264370    .2761319    .4470837          0          1
       wave1 |    264370    .0943261    .2922825          0          1
       wave2 |    264370    .2362863    .4248008          0          1
       wave3 |    264370    .2932027    .4552314          0          1
       wave4 |    264370    .3761849    .4844282          0          1

. desc;

Contains data from C:\matt\publications\ajps4\replication\data_analysis\individual_religious_participation.dta
  obs:       264,370                          
 vars:            50                          10 Oct 2012 07:49
 size:    55,517,700 (86.8% of memory free)
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
              storage  display     value
variable name   type   format      label      variable label
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
year            int    %8.0g       s020       survey year (s020)
ccode           float  %9.0g                  
country         str48  %48s                   country name
attend_religi~s byte   %26.0g      relig      how often do you attend religious services (reverse of f028)
wave            byte   %8.0g       s002       survey wave (s002)
highest_educa~n byte   %8.0g       x025       level of education (x025)
education_rec~e byte   %8.0g       x025r      three levels of education 1 = lower 2 = middle 3 = upper (x025r)
social_class    byte   %18.0g      social     self-reported social_class (reverse of x045)
income_scale    byte   %8.0g       x047       income_scale (x047)
income_level    byte   %8.0g       x047r      three levels of income (x047r)
male            float  %9.0g                  female = 0, male = 1
hdi             float  %8.0g                  human development index
government_re~n float  %8.0g                  government regulation of religion (IRF, ARDA)
social_regula~n float  %8.0g                  social regulation of religion (IRF, ARDA)
government_fa~m byte   %8.0g                  government favoritism of religion (IRF, ARDA)
relfrac_alesina float  %8.0g                  religious fractionalization (alesina)
protestant_al~a float  %8.0g                  percent protestant (alesina)
roman_catholi~a float  %8.0g                  percent roman catholic (alesina)
muslim_alesina  float  %8.0g                  percent muslim (alesina)
taoist_alesina  float  %8.0g                  percent taoist (alesina)
buddhist_ales~a float  %8.0g                  percent buddhist (alesina)
inequality_gini float  %8.0g                  Inequality (gini) by Banones
urbanization    byte   %8.0g                  percentage urban (WDI online)
regulation      byte   %8.0g                  overall regulation of majority religion (RAS)
gir             float  %8.0g                  overall government involvement in religion (RAS)
region          byte   %8.0g                  geographic region (aclp)
polity2         byte   %8.0g                  Polity score
gdp1000         float  %9.0g                  GDP per capita in thousands 1996 PPPs (PWT 6.1)
communist       float  %9.0g                  communist at time of survey
postcommunist   float  %9.0g                  former communist state
protestant_ma~a float  %9.0g                  0 = no, 1= yes
catholic_majo~a float  %9.0g                  0 = no, 1= yes
muslim_majori~a float  %9.0g                  0 = no, 1= yes
lhdi            float  %9.0g                  ln of HDI
lgdp1000        float  %9.0g                  ln of gdp1000
age65           float  %9.0g                  1 = age 65 or over, 0 otherwise
sub_saharan     float  %9.0g                  1 = sub-Saharan, 0 otherwise
south_asia      float  %9.0g                  1 = South Asia, 0 otherwise
east_asia       float  %9.0g                  1 = East Asia, 0 otherwise
se_asia         float  %9.0g                  1 = South-East Asia, 0 otherwise
oceania         float  %9.0g                  1 = Oceania, 0 otherwise
middle_east     float  %9.0g                  1 = Middle East, 0 otherwise
latin_america   float  %9.0g                  1 = Latin America, 0 otherwise
north_america   float  %9.0g                  1 = North America, 0 otherwise
east_europe     float  %9.0g                  1 = East Europe, 0 otherwise
west_europe     float  %9.0g                  1 = West Europe, 0 otherwise
wave1           float  %9.0g                  1 = WVS Wave 1, 0 otherwise
wave2           float  %9.0g                  1 = WVS Wave 2, 0 otherwise
wave3           float  %9.0g                  1 = WVS Wave 3, 0 otherwise
wave4           float  %9.0g                  1 = WVS Wave 4, 0 otherwise
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Sorted by:  

. *     ****************************************************************  *;
. *       Create a panel ID variable.                                     *;
. *     ****************************************************************  *;
. egen idn=concat(year ccode);

. encode idn, gen(id);

. *     ****************************************************************  *;
. *       Get some descriptive information.                               *;
. *     ****************************************************************  *;
. xtset id;
       panel variable:  id (unbalanced)

. xtsum attend_religious_services;

Variable         |      Mean   Std. Dev.       Min        Max |    Observations
-----------------+--------------------------------------------+----------------
attend~s overall |  4.335658   2.558971          1          8 |     N =  249063
         between |             1.292941   1.130832   7.398033 |     n =     180
         within  |             2.226959  -2.062375   11.20483 | T-bar = 1383.68

. xttab  attend_religious_services;

                  Overall             Between            Within
attend_~s |    Freq.  Percent      Freq.  Percent        Percent
----------+-----------------------------------------------------
 never pr |   61125     24.54       177     98.33          25.17
 less tha |   25165     10.10       177     98.33          10.75
 once a y |   19115      7.67       178     98.89           8.17
 other sp |    6957      2.79        83     46.11           6.20
  only on |   36318     14.58       179     99.44          14.90
 once a m |   28333     11.38       180    100.00          11.24
 once a w |   44764     17.97       180    100.00          17.41
 more tha |   27286     10.96       180    100.00          10.29
----------+-----------------------------------------------------
    Total |  249063    100.00      1334    741.11          13.49
                               (n = 180)

. *     ****************************************************************  *;
. *       Replicate Table 1                                               *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. *       Table 1, Model 1                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, vce (cluster id) theta;

Random-effects GLS regression                   Number of obs      =    229262
Group variable: id                              Number of groups   =       164

R-sq:  within  = 0.0000                         Obs per group: min =       409
       between = 0.6675                                        avg =    1397.9
       overall = 0.1757                                        max =      4525

Random effects u_i ~ Gaussian                   Wald chi2(11)      =    571.54
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8581   0.8926     0.9160     0.9451   0.9570

                                   (Std. Err. adjusted for 164 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        lhdi |  -2.994269   .5727989    -5.23   0.000    -4.116934   -1.871604
urbanization |  -.0232234   .0065543    -3.54   0.000    -.0360696   -.0103773
inequality~i |   .0351427   .0077038     4.56   0.000     .0200435    .0502418
government~n |  -.1440336   .0400998    -3.59   0.000    -.2226278   -.0654394
government~m |  -.0086686    .081148    -0.11   0.915    -.1677157    .1503786
social_reg~n |   .0474927   .0318615     1.49   0.136    -.0149546    .1099401
   communist |  -.8958044   .5470372    -1.64   0.102    -1.967978    .1763687
postcommun~t |  -.1577307   .1864004    -0.85   0.397    -.5230687    .2076073
roman_cath~a |   .0139243   .0026817     5.19   0.000     .0086682    .0191804
pr~t_alesina |  -.0009091   .0029357    -0.31   0.757    -.0066629    .0048447
muslim_ale~a |   .0092717   .0039086     2.37   0.018      .001611    .0169325
       _cons |   3.559538   .5784017     6.15   0.000     2.425891    4.693184
-------------+----------------------------------------------------------------
     sigma_u |  .76416729
     sigma_e |  2.2149744
         rho |   .1063652   (fraction of variance due to u_i)
------------------------------------------------------------------------------

.         *     ****************************************************************  *;
. *       Table 1, Model 2                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services income_level lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =    198341
Group variable: id                              Number of groups   =       158

R-sq:  within  = 0.0017                         Obs per group: min =       333
       between = 0.6682                                        avg =    1255.3
       overall = 0.1767                                        max =      4303

Random effects u_i ~ Gaussian                   Wald chi2(12)      =    619.99
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8445   0.8810     0.9118     0.9421   0.9562

                                   (Std. Err. adjusted for 158 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
income_level |  -.1146916   .0195098    -5.88   0.000      -.15293   -.0764531
        lhdi |  -2.955246   .5941836    -4.97   0.000    -4.119824   -1.790667
urbanization |  -.0237365   .0066551    -3.57   0.000    -.0367803   -.0106928
inequality~i |   .0328648   .0077947     4.22   0.000     .0175875    .0481421
government~n |  -.1473613   .0405947    -3.63   0.000    -.2269255   -.0677971
government~m |  -.0313705   .0857119    -0.37   0.714    -.1993628    .1366218
social_reg~n |   .0511581   .0330205     1.55   0.121    -.0135609    .1158771
   communist |  -.8778916   .5398987    -1.63   0.104    -1.936074    .1802904
postcommun~t |   -.095088    .200025    -0.48   0.635    -.4871298    .2969539
roman_cath~a |   .0146465   .0029043     5.04   0.000     .0089542    .0203387
pr~t_alesina |   .0000218   .0031225     0.01   0.994    -.0060982    .0061418
muslim_ale~a |   .0098821   .0039342     2.51   0.012     .0021711     .017593
       _cons |   3.919555    .590088     6.64   0.000     2.763004    5.076106
-------------+----------------------------------------------------------------
     sigma_u |   .7692468
     sigma_e |  2.2101629
         rho |  .10804968   (fraction of variance due to u_i)
------------------------------------------------------------------------------

. *     ****************************************************************  *;
. *       Table 1, Model 3                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services male lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =    229111
Group variable: id                              Number of groups   =       164

R-sq:  within  = 0.0099                         Obs per group: min =       407
       between = 0.6681                                        avg =    1397.0
       overall = 0.1833                                        max =      4525

Random effects u_i ~ Gaussian                   Wald chi2(12)      =    587.28
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8588   0.8931     0.9166     0.9455   0.9572

                                   (Std. Err. adjusted for 164 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
        male |  -.4418885   .0833833    -5.30   0.000    -.6053168   -.2784603
        lhdi |  -3.017599   .5726563    -5.27   0.000    -4.139985   -1.895214
urbanization |  -.0234116   .0065544    -3.57   0.000     -.036258   -.0105651
inequality~i |   .0351345   .0076915     4.57   0.000     .0200593    .0502096
government~n |    -.14332   .0400814    -3.58   0.000     -.221878    -.064762
government~m |  -.0088096   .0813269    -0.11   0.914    -.1682074    .1505882
social_reg~n |   .0477882   .0319219     1.50   0.134    -.0147777    .1103541
   communist |  -.9008678   .5435949    -1.66   0.097    -1.966294    .1645586
postcommun~t |  -.1638218   .1868087    -0.88   0.381    -.5299601    .2023164
roman_cath~a |   .0139954   .0026827     5.22   0.000     .0087374    .0192534
pr~t_alesina |  -.0007085   .0029318    -0.24   0.809    -.0064548    .0050379
muslim_ale~a |   .0093893   .0039218     2.39   0.017     .0017026    .0170759
       _cons |   3.769161   .5768642     6.53   0.000     2.638528    4.899794
-------------+----------------------------------------------------------------
     sigma_u |   .7656807
     sigma_e |  2.2039866
         rho |  .10769395   (fraction of variance due to u_i)
------------------------------------------------------------------------------

.         *     ****************************************************************  *;
. *       Table 1, Model 4                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services age65 lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =    227504
Group variable: id                              Number of groups   =       162

R-sq:  within  = 0.0099                         Obs per group: min =       406
       between = 0.6675                                        avg =    1404.3
       overall = 0.1840                                        max =      4521

Random effects u_i ~ Gaussian                   Wald chi2(12)      =    813.40
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8598   0.9060     0.9173     0.9460   0.9576

                                   (Std. Err. adjusted for 162 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
       age65 |   .7082881   .0464172    15.26   0.000      .617312    .7992642
        lhdi |  -3.112966   .5761392    -5.40   0.000    -4.242178   -1.983754
urbanization |  -.0233895   .0066018    -3.54   0.000    -.0363288   -.0104503
inequality~i |   .0364928   .0077924     4.68   0.000     .0212201    .0517656
government~n |  -.1474298    .040595    -3.63   0.000    -.2269945    -.067865
government~m |  -.0097315   .0811109    -0.12   0.905    -.1687059    .1492429
social_reg~n |   .0483739    .032179     1.50   0.133    -.0146957    .1114436
   communist |  -.8617202   .5450113    -1.58   0.114    -1.929923    .2064822
postcommun~t |  -.1798865    .193289    -0.93   0.352     -.558726     .198953
roman_cath~a |   .0137528   .0028104     4.89   0.000     .0082445     .019261
pr~t_alesina |  -.0012494   .0030832    -0.41   0.685    -.0072923    .0047935
muslim_ale~a |   .0094062   .0039451     2.38   0.017      .001674    .0171384
       _cons |   3.434942   .5827649     5.89   0.000     2.292744     4.57714
-------------+----------------------------------------------------------------
     sigma_u |  .77151307
     sigma_e |  2.2019888
         rho |  .10933759   (fraction of variance due to u_i)
------------------------------------------------------------------------------

. *     ****************************************************************  *;
. *       Table 1, Model 5                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services highest_education lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =    158589
Group variable: id                              Number of groups   =       111

R-sq:  within  = 0.0025                         Obs per group: min =       402
       between = 0.7097                                        avg =    1428.7
       overall = 0.1942                                        max =      4298

Random effects u_i ~ Gaussian                   Wald chi2(12)      =   3726.76
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8458   0.9000     0.9096     0.9423   0.9523

                                   (Std. Err. adjusted for 111 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
highest_ed~n |  -.0512071   .0091439    -5.60   0.000    -.0691287   -.0332855
        lhdi |  -2.571194   .7160954    -3.59   0.000    -3.974715   -1.167673
urbanization |  -.0278836   .0090826    -3.07   0.002    -.0456852   -.0100819
inequality~i |   .0388758   .0075659     5.14   0.000     .0240469    .0537048
government~n |  -.0796584   .0529832    -1.50   0.133    -.1835035    .0241867
government~m |   .0563729   .0946025     0.60   0.551    -.1290446    .2417903
social_reg~n |   .0429852   .0435933     0.99   0.324     -.042456    .1284264
   communist |  -3.330656   .4997971    -6.66   0.000     -4.31024   -2.351072
postcommun~t |  -.1791502   .2552567    -0.70   0.483    -.6794441    .3211438
roman_cath~a |   .0129669   .0034695     3.74   0.000     .0061667    .0197671
pr~t_alesina |   .0006463   .0043567     0.15   0.882    -.0078927    .0091852
muslim_ale~a |   .0056054   .0039635     1.41   0.157    -.0021628    .0133736
       _cons |    3.74891   .8963997     4.18   0.000     1.991999    5.505821
-------------+----------------------------------------------------------------
     sigma_u |  .69759731
     sigma_e |  2.1834624
         rho |   .0926205   (fraction of variance due to u_i)
------------------------------------------------------------------------------

. *     ****************************************************************  *;
. *       Table 1, Model 6                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services social_class lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =    111397
Group variable: id                              Number of groups   =        79

R-sq:  within  = 0.0001                         Obs per group: min =       406
       between = 0.7608                                        avg =    1410.1
       overall = 0.2021                                        max =      3280

Random effects u_i ~ Gaussian                   Wald chi2(12)      =   3856.13
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8468   0.8929     0.9093     0.9421   0.9455

                                    (Std. Err. adjusted for 79 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
social_class |   .0271462    .020854     1.30   0.193    -.0137269    .0680193
        lhdi |  -2.753016   .8659113    -3.18   0.001    -4.450171   -1.055861
urbanization |  -.0234579   .0100846    -2.33   0.020    -.0432233   -.0036924
inequality~i |   .0424689   .0095903     4.43   0.000     .0236723    .0612655
government~n |  -.0737842   .0617949    -1.19   0.232       -.1949    .0473315
government~m |   .0467413   .1126238     0.42   0.678    -.1739973    .2674799
social_reg~n |   .0575284   .0475121     1.21   0.226    -.0355935    .1506503
   communist |  -2.957478    .527407    -5.61   0.000    -3.991176   -1.923779
postcommun~t |  -.0244761   .2910188    -0.08   0.933    -.5948626    .5459103
roman_cath~a |   .0161759   .0032101     5.04   0.000     .0098842    .0224676
pr~t_alesina |    .006433   .0045195     1.42   0.155     -.002425     .015291
muslim_ale~a |    .009722   .0043635     2.23   0.026     .0011697    .0182743
       _cons |   2.579248   .9148268     2.82   0.005     .7862201    4.372275
-------------+----------------------------------------------------------------
     sigma_u |  .69591145
     sigma_e |  2.1734535
         rho |  .09298672   (fraction of variance due to u_i)
------------------------------------------------------------------------------

.         *     ****************************************************************  *;
. *       Table 1, Model 7                                                *;
. *     ****************************************************************  *;
. xtreg attend_religious_services income_level male age65 highest_education social_class lhdi urbanization inequality_gini  
>         government_regulation government_favoritism social_regulation
>         communist postcommunist roman_catholic_alesina protestant_alesina muslim_alesina, cluster(id) i(id) theta;

Random-effects GLS regression                   Number of obs      =     95555
Group variable: id                              Number of groups   =        73

R-sq:  within  = 0.0114                         Obs per group: min =       325
       between = 0.7761                                        avg =    1309.0
       overall = 0.2166                                        max =      2962

Random effects u_i ~ Gaussian                   Wald chi2(16)      =   4847.12
corr(u_i, X)       = 0 (assumed)                Prob > chi2        =    0.0000

------------------- theta --------------------
  min      5%       median        95%      max
0.8167   0.8829     0.8997     0.9342   0.9384

                                    (Std. Err. adjusted for 73 clusters in id)
------------------------------------------------------------------------------
             |               Robust
attend_rel~s |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
income_level |  -.0597257   .0271722    -2.20   0.028    -.1129821   -.0064693
        male |  -.2984654   .1316278    -2.27   0.023    -.5564512   -.0404795
       age65 |   .4632703   .0672574     6.89   0.000     .3314482    .5950924
highest_ed~n |  -.0351441   .0099705    -3.52   0.000    -.0546859   -.0156022
social_class |   .0896118   .0165143     5.43   0.000     .0572444    .1219791
        lhdi |  -2.985337   .8998107    -3.32   0.001    -4.748934    -1.22174
urbanization |  -.0239818   .0099802    -2.40   0.016    -.0435427   -.0044209
inequality~i |   .0387652   .0096241     4.03   0.000     .0199023     .057628
government~n |  -.0899261   .0569186    -1.58   0.114    -.2014845    .0216322
government~m |   .0261181   .1233488     0.21   0.832    -.2156412    .2678773
social_reg~n |   .0755719   .0458068     1.65   0.099    -.0142078    .1653515
   communist |  -2.640071    .545107    -4.84   0.000    -3.708461   -1.571681
postcommun~t |   .2406978   .3531502     0.68   0.496    -.4514638    .9328595
roman_cath~a |   .0201005   .0041443     4.85   0.000     .0119779    .0282232
pr~t_alesina |    .010426   .0053761     1.94   0.052     -.000111     .020963
muslim_ale~a |   .0119488   .0048896     2.44   0.015     .0023654    .0215321
       _cons |   2.662392   .9083363     2.93   0.003     .8820854    4.442698
-------------+----------------------------------------------------------------
     sigma_u |  .64357749
     sigma_e |  2.1631978
         rho |  .08131589   (fraction of variance due to u_i)
------------------------------------------------------------------------------

.         *     ****************************************************************  *;
. *     ****************************************************************  *;
. *       Replication of Table 1 complete                                 *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. *       The End                                                         *;
. *     ****************************************************************  *;
. *     ****************************************************************  *;
. 
end of do-file

. log close
      name:  <unnamed>
       log:  C:\matt\publications\ajps4\replication\data_analysis\individual_religious_participation.log
  log type:  text
 closed on:  11 Oct 2012, 05:42:24
--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
